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Resume 

Soft A un FjTj-ordre maximal d'un corps de quaternions sur F 9 (T) non- 
ramifie a la place oo. Cet article donne un algorithme pour calculer un do- 
maine fondamental de Taction de groupe des unites A* sur l'abre de Bruhat- 
Tits T associee a PGL 2 (F g ((l/T))), Faction qui est un analogue en corps de 
functions de Taction d'un groupe cocompact Fuchsian sur le demi-plan su- 
perieur. L'algorithme donne egalement une presentation explicit du groupe 
A* par generateurs et relations. En outre nous trouvons une borne superieure 
pour le temps de calcul en utilisant que le graphe quotient A*\T est presque 
Ramanujan. 

Abstract 

Let A be a FjTj-maximal order in a division quaternion algebra over F 9 (T) 
which is split at the place oo. The present article gives an algorithm to 
compute a fundamental domain for the action of the group of units A* on 
the Bruhat-Tits tree T associated to PGL 2 (F q ((l/T))). This action is a 
function field analog of the action of a co-compact Fuchsian group on the 
upper half plane. The algorithm also yields an explicit presentation of the 
group A* in terms of generators and relations. Moreover we determine an 
upper bound for its running time using that A*\T is almost Ramanujan. 



1 Introduction 

A major recent theme in explicit arithmetic geometry over Q or more general 
number fields has been the development and implementation of algorithms 
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to compute automorphic forms [Cr[ IDej IGVj IGY[ IStej . More precisely, these 
algorithms compute the Hecke action on spaces of modular forms for a given 
level and weight. Typically these algorithms proceed in three steps: (i) a 
combinatorial or geometric model is provided in which one can compute the 
Hecke action; (ii) on the model one performs some precomputations such 
as the computations of ideal classes of a maximal order in a quaternion 
division algebra, or the computation of a fundamental domain; (iii) on the 
data provided by (ii) one implements the Hecke action. 

The present article is concerned with an analogous algorithm over func- 
tion fields whose ultimate goal is the computation of Drinfeld modular forms 
as well as automorphic forms. For GL 2 over function fields, (i) and (ii) were 
solved in |Telj and \Te2\ IGNj . respectively. Here we will be concerned with 
inner forms of GL 2 that correspond to the unit group of a quaternion division 
algebra split at oo. In this setting an extension of |Telj is part of |Buj . The 
sought-for combinatorial description of the forms to be computed is given in 
terms of harmonic cocycles on the Bruhat-Tits tree which are invariant un- 
der the action of an arithmetic subgroup T defined from the division algebra. 
The main precomputation that makes up step (ii) is that of a fundamental 
domain of T under the action of T. This can be thought of as an analog 
of [Voj . Due to the different underlying geometry the methods employed are 
completely different. 

To describe the output of our algorithm and some consequences note first 
that in our setting of a quaternion division algebra split at oo, the quotient 
r\T is a finite graph. The fundamental domain with an edge pairing that 
we compute consists of the following data: 

1. a finite subtree Y C T whose image Y in T\T is a maximal spanning 
tree, i.e., Y is a tree such that adding any edge of T\T to it will create 
a cycle. 

2. for any edge e of T\T \ Y, an edge e of T connected to Y that maps 
to e and a the gluing datum that connects the loose vertex of this edge 
via the action of T to a vertex of Y. 

Compared to output of |Voj . what we compute is the analog of fundamental 
domain together with a side pairing. As explained in |Sel|. § 1.4], this data 
yields a presentation of the group V in terms of explicit generators and rela- 
tions. Moreover the fundamental domain data computed provides an efficient 
reduction algorithm on the tree: to any edge it computes its representative 
in Y', by which we mean the union of Y with the edges in fl2]). Reinterpreted 
in terms of group theory, a fundamental domain with an edge pairing yields 
an efficient algorithm to solve the word problem for V. 
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Observing that a finite cover of T\T is a Ramanujan graph, yields a 
bound on the diameter of T\T. This in turn we use to bound the complexity 
of our algorithm, to bound the size of Y', and to bound the size of the 
representatives of Y' in terms of a natural height on the 2 x 2-matrices 
over the function field completed at oo. Our main result is therefore the 
existence of an effective algorithm together with precise complexity bounds. 
An implementation can be obtained on request from the second author. 

A theoretical result on the size of a minimal generating set for T and the 
(logarithmic) height of its generators was obtained by different methods in 
|Pa2] . The bound on the height here is better by a factor of 2. The results 
in |Pa2j also prove the existence of an algorithm to compute a fundamental 
domain. An implementation or a detailed analysis of that algorithm have 
not yet been carried out. Both |Pa2j and the present article rely heavily on 
|Pal] by Papikian. 

We conclude by a short overview of the article: in Sections HJ [3] and HI we 
recall basic notions and results from graph theory, on the Bruhat-Tits tree 
and from the theory of quaternion algebras. Section [5] introduces the main 
object of this article, the action of T on the Bruhat-Tits tree T, and states 
relevant results on the resulting quaternion graph T\T. Our basic algorithm 
is presented in Section except that we use one unproved subroutine which 
is the content of the following Section [7J The final Sections [S] and M present 
on the one hand the applications of the algorithm to the presentation of T in 
terms of generators and relations and to the word problem, and on the other 
hand the complexity analysis of the algorithm based on the fact that T\T 
has a finite cover that is Ramanujan. 

Acknowledgments: For several useful discussions we wish to thank 
Mihran Papikian and John Voight. We also want to heartily thank the 
anonymous referee whose many comments and suggestions greatly improved 
the readability of the present article. During this work, the authors were sup- 
ported by the Sonderforschungsbereich/Transregio 45 Periods, Moduli Spaces 
and Arithmetic of Algebraic Varieties and by the DFG priority project SPP 
1489. The implementation of the algorithm is based on the computer algebra 
system Magma |BCPj . 

Notation 

Throughout this article K = F g (T) will denote the rational function field 
over F g . As usual, the infinite valuation on K is defined by Voo(-) = 
deg(g) - deg(f) for f,g e A = ¥ g [T],g ^ and Voo (0) = oo. Then vr = l/T 
is a uniformizer for v^, the corresponding completion of K is = F 9 ((7r)) 
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and we write for its ring of integers. 

Remark 1.1. The restriction that q be odd is for the sake of simplicity of 
exposition. To treat the case that q is even, one needs to redo all of Section 0] 
with very little overlap with the odd case. Some changes are also necessary 
in Section [7J All other results hold independently of q being even or odd. 
We have implemented our algorithm also for even q. Theoretical discussions 
on q even are sketched in [BuJ. 

2 Notions from graph theory 

Definition 2.1. 1. A (directed multi-) graph Q is a pair (V((/), E(C/)) where 
V(£?) is a (possibly infinite) set and E(£?) is a subset of V(£7) x V(£7) x 
Z> such that 

(a) if e = (v, v', i) lies in E(C/), then so does its opposite e* = (v 1 , v, i), 

(b) for any (v,v') G V(£) x V(£), the set {i G Z> | (v,v',i) G E((?)} 
is a finite initial segment of Z>o of cardinality denoted by n v y , 

(c) for any v G V(£), the set Nbs(w) := {V G V(0) | {v, v', 0) G E((?)} 
is finite. 

An element v G V(£?) is called a vertex, an element e G E(C?) is called an 
(oriented) edge and an element in V(Q) U E((/) is called a simplex. For each 
edge e = (v,v',i) G E(C?) we call o(e) := f the origin of e and t(e) := v' the 
target of e. If there is only one edge between vertices v, v' of Q we simply 
write (v,v') instead of (v,v',0). Two vertices v,v' are called adjacent, if 
{t>,i/} = {o(e),t(e)} for some edge e. 

Definition 2.2. An edge e with o(e) = t(e) is called a loop. A vertex t> with 
# Nbs(f ) = 1 is called terminal. 

Let f , f ' G V(C?). A pai/i from t> to v' is a finite subset {ex, . . . , e^} of E(C?) 
such that t(ej) = o(ej + i) for all i = 1, . . . , k — 1 and o(ei) = v,t(ek) = v'. 
The integer is called the length of the path {ei, . . . , e^}. The distance from 
f to f', denoted d(v,v'), is the minimal length among all paths from v to 
v' (or oo if no such path exists). A path {ei, . . . , e^} from v to v' without 
backtracking, i.e., such that for no i we have = e*_ 1; is called a geodesic. 
Note that the length of a geodesic need not be d{v,v') but that d{v,v') is 
attained for a geodesic. 

A graph Q is connected if for any two vertices f , t> ' G V(C?) there is a path 
from f to v'. A cyc/e of ^ is a geodesic from some vertex v to itself. Therefore 
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a loop is a cycle of length one. A graph Q is cycle-free if it contains no cycles. 
A tree is a connected, cycle-free graph. If Q is a tree, then any two vertices 
of Q are connected by a unique geodesic. 

A subgraph Q' C is a graph £?' such that V(£') C V(£) and E(£') C 
E(£?). Any subgraph 5 C £ which is a tree is called subtree. A maximal 
subtree is a subtree which is maximal under inclusion among all subtrees 

of g. 

The degree of v G V(£7) is 

deg(w) := deg g (w) := ra„ it(e ). 

eeE(S):o(e)=i> 

Thus f is terminal precisely if deg(t>) = 1. A graph Q is called k-regular if 
for all vertices v G V(C?) we have deg(t> ) = k. 

A graph Q is finite, if # V(£?) < oo. Then also # E(£?) < oo, since deg(v) 
is finite for all t> G V(G)- The diameter of a (finite) graph £/ is 

diam(^) := max d(v,v'). 

v,v'EV(G) 

Definition 2.3. The first Betti number h\{Q) of a finite connected graph is 

h(G) :=£?SQ-#V(0) + 1. 

A graph defines an abstract simplicial set. Its geometric realization is 
a topological space \Q\. For finite graphs one has hi(Q) = diniQ Hi(\Q\, Q), 
i.e., the Betti number counts the number of independent cycles of Q. 



3 The Bruhat-Tits tree 

In this section, we recall the definition of the Bruhat-Tits tree for the group 
PGL^i^oo)- It is an important combinatorial object for the arithmetic of K. 
The material can be found in |Selj . 

One defines a graph (V(T),E(T)) as follows: Two Coo-lattices L, L' C 
are called equivalent if there is a A G with V = XL. The set V(T) is 
the set of equivalence classes [L] of such lattices. The set E(T) is the set of 
triples ([L], [L']) such that L, V C are Coo-lattices with ttL C V C L, or 
equivalently such that V C- L and L/L' = F g as 0ocr m odules. In particular 
there is at most one edge between any two vertices. 

By [Sell § II. 1] the graph T = (V(T), E(T)) is a g + 1-regular tree - recall 
that q is the cardinality of the residue field of K^. The group PGL^-fToo) 
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acts naturally on lattice classes by left multiplication (g, [L]) i— > [gL]. This 
induces an action on T. Because the definitions are a special case of a general 
construction, one calls T the Bruhat-Tits tree for PGL2(-ft' o)- 

Let ei = (1,0)* and e 2 = (0, 1)* be the standard basis of K^, thought 
of as column vectors. Write for O^ei © Oooe 2 . The following result is 
well-known and straightforward, using the class equation and the transitive 
action of GL 2 (K OQ ) on bases of K^. 

Proposition 3.1. The map 

<p : GL 2 {K O0 )/GL 2 {0 O0 )Kl -> V(T) 

A H- [AOl] 

is a bijection of left GL^K^-sets. 

The map <fi of the above proposition allows us to represent vertices of T 
by elements of GL^-Koo)- Row-reduction to the echelon form of a matrix 
yields a standard representative in GL^fToo) as expressed by the following 
result. 

Lemma 3.2. Every class of ' Gh^K^) / 'GL 2 (Ooo)^m has a unique represen- 
tative of the form 



with nfZ and g G i^ 0O /7r n C o; called its vertex normal form. 

We also need a criterion for adjacency for matrices in vertex normal form. 
Lemma 3.3. Consider the two matrices in vertex normal form 

A:= (*o ?)' b: =ct 9+ r) 

with n G Z,a G ¥ q , g G K 00 /7r n C 00 and Ze£ Li and L 2 fre the two lattices 

L x :=AOl,L 2 :=BOl. 
Then L\ D L 2 and L\jL 2 = F 9 . 

Remark 3.4. Lemma 13.31 only displays q vertices adjacent to [L x ]. The 

Ar" -1 q\ 

missing one is the class of I y J O^. with g now being replaced its 

class in K OQ /-n n ~ 1 O oa . 
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with q G F* , 0i G F g 



Figure 1: The tree T with the corresponding matrices 

Figure [TJ below illustrates the tree together with the matrices in normal 
form corresponding to vertices. The identification is clear from the previous 
lemma. Note that each line in the picture symbolizes a whole fan expanding 
to the right. The elements a G F* /3 G ¥ q agree on each fan. 

fn n \ 

Write L(n,g) for the Coo-lattice (t>i, ^2)0^ where Vl = \ q J anc ^ 1,2 = 

. Note that L(n, g) = L(n, g') if and only if g = g' (mod ^Ooo). 

Remark 3.5. For n G Z, g G we define 

5 := deg n (^) := min{z G Z>o | g G Tt^O^}. 
Then the path from L(n,g) to L(0, 0) in T is 

L(n,g) - - L(n- l,g) - - - L(n-S,g) = L(n-6,0) - 

- L(sgn(n-5)-(|n-5|-l),0) - - L(sgn(n-5) -1, 0) - - L(0, 0) 

In particular the distance between L(n,g) and L(0, 0) is deg n (g) + |n — 
deg n (»|. 

4 Quaternion algebras 

We recall standard facts on quaternion algebras over K = F 9 (T) and over 
completions of K, and on orders over A = F 9 [T]. We assume throughout 
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that q is odd. Our basic references are p21 Kap. IX] and [VI]. Many results 
stated are true more generally. However, we confine ourselves to the case at 
hand. 

A quaternion algebra over a field F is a central simple algebra of dimen- 
sion 4 over F. It is either isomorphic to M 2 (F) or a division algebra. One 
has the following well-known construction of quaternion algebras. 

Construction 4.1. For a,b G F* one defines (j=r) as the K-algebra with 
F -basis l,i,j,ij and relations i 2 = a,j 2 = b,ij = —ji. 

The relations can be expanded to a 4 x 4 multiplication table for the given 
F-basis of (^r). One shows that (^) defines a quaternion algebra over F, 
and that conversely any quaternion algebra over F can be obtained via this 
construction for a suitable choice of a, b G F*. 

Among other things, a quaternion algebra D over F carries a reduced 
norm map nrd: D — >■ F which defines a quadratic form on D. For D = (%■') 
the reduced norm has the explicit expression 

nrd(7) = 7-7 = 7? - a 2 7 2 - & 2 7^ + a&74 

for any 7 = 7l + 72 i + 73 j + 74 r/ g -D with (71,72,73,74) G F 4 . If D is 
isomorphic to M 2 (F), then nrd is simple the determinant map det : M 2 (F) — > 
F. 

Let now D denote a quaternion algebra over K. Then D p := D ®k Kp is 
a quaternion algebra over the completion K p for any place p of K. 

Definition 4.2. D is ramified at p if and only if D v is a division algebra. 

Definition 4.3. The Hilbert symbol of a pair (a, b) G F 2 at a place p is 



(a,b) 



-1 f^r) is unramified at p 
1 (^) is ramified at p. 



Definition 4.4. For a G A and w an irreducible element of A, the Legendre 
symbol of a at w is 



1 a zuA and a is a square modulo w 
-1 a is a non- square modulo w 
a G wA. 



By adaptating to the function-field situation the proof of [Se2l Ch. Ill, 
Thm. 1], the following result is straightforward. 
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Proposition 4.5. Suppose q is odd. Write p = (to) and let a = oo a u,b = 
zo^v with u, v G 0* K , a, (5 G Z and let e(p) := degfw) (mod 2). TTien 




Let i? denote the set of all ramified places of D. Then [VT1 Lem. III. 3.1 
and Thme. III. 3.1] yields the following. 

Proposition 4.6. The cardinality of R is finite and even and D is up to 
isomorphism uniquely determined by R. The set R is empty if and only if 
D = M 2 {K). 

The ideal r := ripeR p °f ^ ^ s called the discriminant of D. Let r G A be 
the monic generator of r. 

Assumption 4.7. For the remainder of this article, we assume that D is 
a division quaternion algebra which is unramified at oo, , i.e., that D is an 
indefinite quaternion algebra over A. We also fix an isomorphism = 
M 2 (K oc ). ' 

Let A be an order of D over A. It is free over A of rank 4 and so we may 
choose a basis fi, . . . , f±. The ideal generated by 

disc(/i, . . . , f^ := det(trd(/i/_ 7 )) iJ= i,... i 4 

is independent of the chosen basis. By [VU Lem. 1.4.7], this ideal is a square 
and one defines the reduced discriminant rdisc(A) of A as the square root 
of this ideal. One deduces a criterion for an order to be maximal, see [VU 
Cor. III.5.3]. 

Proposition 4.8. An A-order A is maximal in D if and only z/rdisc(A) = r. 

Since D is split at infinity and K has class number 1, [Vl] Cor. III. 5. 7] 
yields: 

Proposition 4.9. All maximal A-orders A in D are conjugate under D* . 

Let T := A* be the group of units of of a maximal order A. By what 
we have said so far, Y depends uniquely up to conjugation on D, i.e., on K 
and R. 

From global to local compatibilities and explicit local results, one deduces 
the following assertions. 

Lemma 4.10. 1. The reduced norm nrd maps A to A. 
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2. r = {7 G A I nrd( 7 ) G F*}. 

3. The embedding 1: D ^ = M^K^) restricts to a group monomor- 



The following result is well-known. In lack of an explicit reference, we 
shall give a proof. 

Proposition 4.11. Via 1 the group T is a discrete subgroup of GL 2 (K 00 ) . 

Proof. The open sets {1 + ir n M 2 (0 O0 ) \ n G N} form a basis of open neigh- 
borhoods of 1 in GL 2 (K 00 ). After shifting by 1 it suffices to show that 
A n M 2 {0 00 ) is finite, or in other words that A is discrete in M 2 (K oc ). 

To see the discreteness, let T> be the unique locally free coherent sheaf of 
rings of rank 4 over Pp such that A = T(Ap , T>) and such that the completed 
stalk at infinity satisfies = M 2 (0 00 ). Then A n M 2 (0 00 ) = H°(P$ ,V). 
By the Riemann-Roch Theorem, this is a finite-dimensional F q - vector space. 

Alternatively, for D and A constructed later in Propositions 14.151 and 
I4.16[ and the embedding from Lemma I4.17[ the discreteness can be verified 
explicitly, by proving that an A-Basis of A maps to a i^oo-basis of D^. □ 

Given an even set R of finite places of K at which D is ramified, the 
algorithm described in Sections M and [7] will be based on a concrete model 
for (D, A). In the remainder of this section we describe such a model. It will 
consist of an explicit pair (a, b) G K* such that D = {^) and an explicit 
basis of a maximal A-order A of (%r) . 

Let / > 2 be even and let R be a set of Z-distinct prime ideals {pi, . . . , p;} 
of A. Denote by Wi the unique monic (irreducible) generator of pj. Set 
r := Yli w i an d t := Yli Pi where the index i ranges over 1, . . . , / . 

Lemma 4.12. There is an irreducible monic polynomial a G A of even degree 
such that 



phism 





(4.1) 



Any such a also satisfies (^) = 1. 



Proof. Choose any a G A such that 
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for all i. This can be done using the Chinese remainder theorem. By the 
strong form of the function field analogue of Dirichlet's theorem on primes 
in arithmetic progression, \Ko\ Thm. 4.8], the set {a + rb \ b G A} contains 
an irreducible monic polynomial a of even degree. Since a = a (mod wi) we 
have 

for all i. By quadratic reciprocity, |Rot Thm. 3.3], we deduce 

= (_ 1 )^deg a de groi ^O^j = _ x 

since deg(a) is even. But then because I is even, we find 



nei =(-d'=i. 

1=1 



□ 



Remark 4.13. In practice a is rapidly found by the following simple search: 
Step 1: Start with m — 2. 

Step 2: Check for all monic irreducible a E A of degree m whether C^~\ 
I for all 1 < % < I. 

Step 3: If we found an a then stop. Else increase m by 2 and go back to 
Step 2. 



In the function field setting |MSj gives an unconditional effective version 
of the Cebotarov density theorem. This allows us to make Lemma 14.121 
effective, i.e., to give explicit bounds on deg(a) in terms of deg(r). 

Proposition 4.14. Abbreviate d := deg(r). The following table gives upper 
bounds on d a '■= deg(a) depending on q and I: 





q = 3 


Q = 


5,7 


Q = 


= 9 


q>ll 




1<A 


1 = 6 


8 < / 


I < 6 


8 < / 


I < 4 


6 < I 


1 = 2 


4 < I 


d a < 


d + 7 


d + 5 


d+1 


d + 3 


d+1 


d + 3 


d + 1 


d + 3 


d+1 



A basic reference for the results on function fields used in the following 
proof is [StiJ. 
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Proof. Let K' := K{y/W[, . . . , y/wi). Then K' /K is a Galois extension with 
Galois group isomorphic to {±1}' with {±1} = Z/(2); it is branch locus in 
K is the divisor D consisting of the sum of the (wk) and (possibly) oo; the 
constant field of K' is again ¥ q . Denote by g' the genus of K' and by D' the 
ramification divisor of K' / K. The ramification degree at all places is 1 or 2 
and hence tame because q is odd. It follows that deg(D') = #G/2 ■ deg(D). 

Let ir(k) denote the places of K of degree k; let nc{k) denote the places 
p of K of degree k for which Frob p = (—1, . . . , —1) G {±1} Z . Note that the 
elements of irc{k) are in bijection to the monic irreducible polynomials a of 
degree k which satisfy the conditions (14.ip . The following two inequalities 
are from [MS[ Thm. 1 and (1.1)] and the Hurwitz formula, respectively: 

i i k/2 k/2 i 

km - | < 2<j w V + 2 V + (1 + *) d ^ D '>- (4 ' 2 > 

k /2 

|g fc + l-A;7r(A;)| <2g'^— . (4.3) 

2^ = -2#G + deg(D') + 2. (4.4) 
After some manipulations one obtains 

1 qk + 1 g fc/ % deg(^) 2 1 deg(L>) 

To ensure that the right hand side is positive for some (even) k, it thus 
suffices that 

f(k) :=g fe -g fc/2 (2^ 1 (deg(r) + 5) + 2)-(A; + l)2 2i - 1 (deg(r) + l) > 0. (4.5) 

We know that / is the number of prime factors of r and hence that I < deg(r). 
There are at most q places of degree 1 and so for small q such as 3, 5, 7, already 
for small I the degree of r must be quite a bit larger than I. For instance 
if I > 7 and q = 3, then deg(r) > 3/ — 9. Using these considerations and 
simple analysis on f(k), it is simple if tedious to obtain the lower bounds in 
the table. We leave details to the reader. □ 



Proposition 4.15. For a as in Lemma the quaternion algebra D := 

(tk") ^ s ram 'i'fi e d exactly at R. 

Proof. The proof is an immediate consequence of Proposition 14.51 □ 

Since r is a square modulo a, there are e, v G A with deg(e) < deg(a) 
and e 2 = r + vol. 
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Proposition 4.16. A := A + Ai + Aj + A^p- is a maximal A-order of D. 

Proof. This follows easily from Proposition 14.81 by computing the discrimi- 
nant of the A-basis given. □ 

Since a has even degree and is monic, there exists a square root of a in 
Koo. We choose one and denote it by y/a. The following result provides an 
explicit realization for the embedding in Lemma [4.10( c). 

Lemma 4.17. The K -algebra homomorphism i : D — )■ M^(K oo) defined by 

fa \ , . fO 1 
„ and j i y I _ 

—\/ot) r 



i i — y I \^ ^ j— J and j 4 (" ^ 1 induces an isomorphism D <&k K t 



oo — 



M 2 (K OQ ). 

Proof. One verifies t(i) 2 = a, i(j) 2 = r and = —L{j)i{i) by an explicit 

calculation. □ 



5 Facts about quaternion quotient graphs 

In Section [3] we have described the natural action of GL 2 (-ft'oo) o n the Bruhat- 
Tits tree T. In the previous section, starting from D as in Assumption 14.71 
we have produced a discrete subgroup Y C GL 2 (i^oo), the unit group of a 
maximal order. In this section we gather some known results about the in- 
duced action of V on T and the quotient graph T\T. We mainly follow |Pal] . 

Lemma 5.1 ( {Sell Cor. to Prop II. 1]). Forv G V(T) and 7 £ T, the distance 
d(v, jv) is even. 

Proposition 5.2 ( |Pal| Lem. 5.1]). The graph T\T is finite graph. 

Proposition 5.3 f jPaTl Prop. 5.2]). Let v G V(T) and e G E(T). T/ien 
T„ := Stabr(f ) is either isomorphic to ¥* or ¥* 2 . T e := Stabr(e) is isomor- 
phic to ¥*. 

Note that the scalar matrices with diagonal in F* are precisely the scalar 
matrices in T. Clearly they act trivially on T. Hence a stabilizer of a simplex 
is isomorphic to F* precisely if it is the set of scalar matrices with diagonal 
in F*. 

We define T to be the image of T in PGL2(-K'oo) - after what we have just 
seen we have T = T/¥*. Then T v := Stabf(f) is either trivial or isomorphic 
to F* 2 /F^ = Z/(q + 1) and f e := Stab f (e) is always trivial. 

Definition 5.4 ( [Sel| II. 2. 9]). We call a simplex t projectively stable if T t is 
trivial and projectively unstable otherwise. 
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Corollary 5.5. Let v G V(T) be projectively unstable. Then T v acts transi- 
tively on the vertices adjacent to v. 



Let 

if some place in R has even degree, 

1 otherwise 



odd(R) :-- 
and let 



g(R) := 1 + (n(* - 1)) - -^ T 2^- 1 odd(i?) 

where q p = g deg ( p ). Let n : T — >■ T\T be the natural projection. 
Theorem 5.6 ([PaTJ Thm. 5.5]). 1. The graph T\T~ has no loops. 

2. h l (T\T)=g(R). 

3. For v G r\T and v G vt" 1 (w) we /iave: 

(a) v is projectively stable if and only if v has degree q + 1. 

(b) v is projectively unstable if and only if v is terminal. 

4- Let V\ (resp. V q+ \) be the number of terminal (resp. degree q + 1 ) 
vertices ofY\T. Then 

y = 2# R " 1 odd(R) and V q+l = ——(2g(R) - 2 + VA. 

q — 1 

6 An algorithm to compute a fundamental 
domain 

Let the notation T, T be as in the previous section. 

Definition 6.1 ([SeTJ § 1.3]). Let G be a group acting on a graph X . A tree 
of representatives of X (mod G) is a subtree S <Z Q whose image in G\X is 
a maximal subtree. 

The following definition is basically |Sel[ § 1.4.1, Lem. 4], see also [Sel, 
§ 1.5.4, Thm. 13]. Note that (a) differs from [Ml § 1.4.1, Def. 7]. 

Definition 6.2. Let G be a group acting on a tree X. 

1. A fundamental domain for X under G is a pair (S,y) of subgraphs 
S C y C X such that 
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(a) S is a tree of representatives of X (mod G), 

(b) the projection E(y) — > E(G\X) is a bijection, and 

(c) any edge of y has at least one of its vertices in S. 

2. An edge pairing for a fundamental domain y of X under G is a map 

PE := PE (5 , y) := {e G E(3>) \ E(«S) | o(e) G 5} -> G : e ^ p e 

such that g e t(e) G V(<S). We write PE for paired edges. To avoid 
cumbersome notation, we usually abbreviate PEy s by PE. 

3. An enhanced fundamental domain for X under G consists of a funda- 
mental domain, an edge pairing and simplex labels G t '■— Stabe(t) for 
all simplices t of 3^- 

An edge pairing encodes that under the G-action any e = (v, v') G PE 
is identified (paired) with ge = (g e v,g e v') when passing from X to G\X. 
Because X is a tree and the image of S in G\X is a maximal subtree, each 
edge in E(^) \ E(«S) has exactly one of its vertices in V(«S) and therefore 
PE contains exactly those edges of E(^) \ E(<S) pointing away from S. An 
enhanced fundamental domain is a graph of groups in the sense of [Sell 1.4.4, 
Def. 8] realized inside X. Given a fundamental domain with an edge pairing 
the tree S can be recovered from y and PE. 

Remark 6.3. If one barycentrically subdivides T, an alternative way to 
think of an edge pairing is that it pairs the two half sides [o(e), \o(e) + \t(e)] 
and [g e t{e),g e {\o{e) + \t{e))\. 

It will be convenient to introduce the following notation: 

Definition 6.4. For any group G acting on a set X we define a category 
Cg(X) whose objects are the elements of X and whose morphism sets are 
defined as 

Hom G (x,y) := {7 G G | gx = y} C G. 
for x, y G X. The composition of morphisms is given by multiplication in G. 

In particular Endc(a;) := Hom G (x,x) = Stabc(x). 

For the remainder of this section, we assume that Homr(f,w) can be 
computed effectively for all v,w G V(T). This will be verified in Section [71 

Algorithm 6.5. (Computation of the quotient graph) 

Input: A subgroup F C GLi2(K 00 ) for which there exists a routine for com- 
puting Honir(t', v') for all v, v' G V(T) which are equidistant from [L(0, 0)]. 
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Output: A directed multigraph Q with a label attached to each simplex. The 
label values on edges are either (e, 1) (preset), or (e, —1), or a pair (e, g) with 
e G E(T), g G r. The label values on vertices are either (v, 1) (preset) or 
(v, G) for v G V(T) and G dY a finite subgroup. 

Algorithm: 

1. Set v = [£(0,0)]. If #Endr(w ) = q 2 ~ 1, replace v by [£(1,0)]. If 
after replacement we still have ^Endr(fo) = q 2 — 1, then terminate 
the algorithm with the output the connected graph on 2 vertices and 
one edge and with vertex labels Endr(f) for each of the two vertices v. 

2. Initialize a graph Q with V(G) = {v } and E(£?) = 0. Also, initialize 
lists L : = (e G E(T) | o(e) = t> ), the edges adjacent to v , and L' := 0. 
All vertices f of T are given by a matrix in vertex normal form vnf(v). 

3. While L is not empty: 

(a) For % = 1 to #L do: 

i. Let e = (i>,i>') be the ith element in L. 

ii. Compute Endr(f')- 

iii. If #End r «) = Q 2 - 1 then: 

A. Add the vertex v' to V(^) and e and e* to E(^). 

B. Store (v ', Endr(f ')) as a vertex label for v' . 

C. Remove e from L. 

iv. If # Endr(f') — q — 1, then for all j < % do the following: 

A. Let e' = (w,w') be the jth element in L. 

B. Compute Hom r (f', w'). 

C. If Hom r (tj',«j') 7^ 0, then do the following 

• Add an edge e' from v to w' to E(^), as well as its 
opposite. 

• Give e' the label (e,g e ) for some g e G Homr(f',w') and 
give e'* the label (e, —1). 

• Remove (f,f') from L and set j := i. 

• Remove («/, vnf {g e v)) from L'. 

• If now degg(-uZ) = q + 1, then remove («;, it;') from L. 

D. Continue with the next j. 

v. If at the end of the j-loop we have j = i, then: 

A. Add v' to V(£7), add e and e* to E(^). 

B. For all adjacent vertices w ^ v of v ' in T add (f',w) to 
L'. 
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(b) Set L := I! and V := 0. 



4. If L is empty, return Q. 

Remark 6.6. One could randomly choose a vertex [L(n, g)} as Vq and replace 
it by [L(n + 1, g)], if it is projectively unstable. In this case, one would need 
to change the input of Algorithm 16.51 accordingly. 

Remark 6.7. The vertex label (v, 1) is used at all projectively stable vertices. 
For these, the stabilizer is the center of GL 2 (i^oo) intersected with T. There 
is no need to store this group each time. The same remark applies to all edge 
labels (e, 1). 

A maximal subtree S of Q consists of all vertices and those edges with 
edge label (e, 1). It is completely realized within T ■ 

The edges with label (e, g) are the edges which occur (ultimately) in PE. 
The edge label (e, — 1) indicates that the opposite edge has a label (e, g). 
It is clear that the vertex and edge label allow one to easily construct an 
enhanced fundamental domain (S, y) with an edge pairing and labels for the 
action of T on T. 

Theorem 6.8. Suppose T from Algorithm \6.5\ satisfies the following condi- 
tions: 

1. d(v, gv) is even for all g e Y , v e V(T), 

2. for simplices t of T either T t is trivial, or t is a vertex and T t = 

z/Gz + i), 

3. T\T is finite. 

Then Algorithm \6.5\ terminates and computes an enhanced fundamental do- 
main for T under T. 

By the results in Section^ hypotheses (a) -(c) are satisfied ifY is the unit 
group of a maximal order of a quaternion algebra D as in Assumption ^ . 7[ 

Remark 6.9. Let us comment on the hypotheses made in Theorem 16.81 so 
that Algorithm l6.5l terminates: Following the example set by the number field 
case, it seems natural to consider the following situation: Let C be a smooth 
projective geometrically connected curve over ¥ q , let S be a finite set of closed 
points of C, set A = T(C \ S, Oq) and let K = Q(A) be the fraction field of 
A. Let furthermore D denote a division algebra over K which is ramified at 
all but one point oo of S and let A be an A-order. As can be deduced in this 
situation from [Vi] by an argument similar to Proposition 14.111 the group 
A* modulo its center acts discretely on the Bruhat-Tits tree for PGL 2 (-ft'oo)- 
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We expect but have not checked that hypotheses (l)-(3) of Theorem 16.81 are 
always met in this situation. What is missing in this general situation is an 
explicit algorithm to compute Hom\*(v,v'). For this, see Remark 17.51 

Proof of Theorem \6.8i . Let Q be the output of Algorithm 16.51 We show that 
any two distinct simplices of Q have labels (t, ?) and (f, ?) with i! Ft and 
that for all simplices t of T there is a simplex of Q whose label is (f, ?) for 
some t' G Ft. 

For the first assertion, let Vi,v 2 G V(T) be distinct first entries in labels 
of vertices of Q and suppose that jvi = v 2 for some 7 G T \ T Vl . We seek a 
contradiction. In a first reduction step we show that we may assume that v\ 
is projectively stable: So suppose v\ is projectively unstable. Then since 

Stabr(f2) = 7 Stabr(f i)7 _1 , (6.1) 

also v 2 has to be projectively unstable. Hence both v\ and v 2 are terminal 
vertices in Q. Let v[ and v' 2 be their unique adjacent vertices in Q. Since v[ 
is adjacent to v%, it follows that , yv' l is adjacent to 7-17 = v 2 . By condition 
(b) the stabilizer Stabr(t>2) acts transitively on the vertices adjacent to v 2 . 
Hence there exists 7' G Stabr(t>2) with 

iiv'x = v 2 , (6.2) 

and so v[ and v' 2 are T-equivalent. If v[ and v' 2 were also projectively un- 
stable and therefore terminal vertices in Q, then, since Q is connected, Q 
would have to be the graph consisting of the two vertices vi, v 2 and one edge 
connecting them. This contradicts condition (a). Therefore v[ and v' 2 must 
be projectively stable and T-equivalent. To conclude the reduction, observe 
that we cannot have v[ = v 2 , since in this case we must have 7^ G ¥* from 
(16.21) . But 7 r 7 maps v\ to v 2 and this would contradict V\ 7^ v 2 . 

Now suppose v\ is projectively stable. Then by equation (16.11) so is v 2 . Let 
v be the initial vertex of the algorithm and let i\ = d(v, v±) and i 2 = d(v, v 2 ). 
We prove the assertion by induction over i\\ If i\ — 1 then also i 2 = 1 because 
of condition (a). Hence the vertices V\ and v 2 both have the same distance 
1 from v and since Hornet"!, v 2) = q — 1, Algorithm 16.51 with the first choice 
of L rules out that they both lie in Q. This is a contradiction. The same 
reasoning rules out i\ = i 2 for any i±,i 2 > 1. 

Suppose %x > 1. By condition (a) and the previous line we may assume 
ii = i 2 + 2m for some m G Z>i. Let v[ be the vertex on the geodesic from V\ 
to v so that d{y , v[) = i\ — 1. Then by the construction of Q we have v [ G Q. 
The vertex , jv[ is adjacent to 717 = v 2 . Now observe that 7^ does not belong 
to Q because otherwise we could apply the induction hypothesis to v[,-yv[, 
using d(v[, v) = ii — 1 and 0^(7^, v) < i 2 + 1 to obtain a contradiction. 
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It follows that v' 2 := 71^ Since by construction the geodesic from v 
to t>2 lies on ^ we have c?(f 2 , u) = i% + 1- Now by the algorithm that defines 
Q the vertex v' 2 must be equivalent to a vertex of distance — 1, i.e., there 
are 7' G r, v 2 ' G with d(v 2 ,v) = 12 — 1 such that t> 2 ' = 7't> 2 . But then we 
apply the induction hypothesis to v' x ,v 2 and again obtain a contradiction. 
This concludes the proof of the first assertion for vertices. 

Suppose now that e = (v , Vi), e'(v' , v'^ occur as first entries in E((?), lie 
in the same T-orbit, are distinct and occur within edge labels of Q. Let 7 
be in V with e' = 7c Note that not all the vertices V{ and v\ must occur in 
vertex labels from Q but each edge must at least have one vertex that does 
- see step (c)(i)4.C. Suppose after possibly changing the orientation of edges 
and the indices that vq has minimal distance from v. By construction of Q 
the vertex Vq occurs in a vertex label. If v' = 71*0 occurs in a vertex label 
of Q, then by the case already treated, we must have t> = v' Q . Since e^e' 
it follows that v$ is projectively unstable. But then the algorithm does not 
yield an edge starting at v and ending at a vertex v\ with d(v, Vx) > d(v, v ). 
This is a contradiction. 

It follows that v' = 7t>o does not occur in a vertex label. Hence v[ 
must occur in a vertex label. By essentially the argument just given, v\ can 
also not occur in a vertex label. Hence (e, 7) must be an edge label and 
moreover d(v,v[) = d(v,vo) + 1 = d(v,v' ) — 1. But then in step (c)(i)4.C of 
Algorithm 16.51 the edge e' must have been removed from the list V and so it 
cannot occur in a label of an edge of Q. 

We finally come to the second assertion: By construction, Q defines a 
connected subgraph of T\T, since we already showed that there are no T- 
equivalent simplices in Q. Moreover, at any vertex of this subgraph the 
degree within Q and within T\T is the same. Hence Q defines a connected 
component of Y\T. But T and hence T\T are connected and thus Q = 

r\r. □ 

We further describe an algorithm to compute for any v' G V(T) a T- 
equivalent vertex v" G Q. This can be done in time linear to the distance 
from v' to Q. For this algorithm we need the stabilizers of the terminal 
vertices of Q and the elements 7 G Homr(ui, Vj), which we both stored as 
vertex and edge labels during the computation of Q . We call this algorithm 
the reduction algorithm. We need to be able to do the following: 

1. Find the geodesic from v' to v. This was discussed in Remark 13.51 

2. Determine the extremities of a given geodesic in Q. Since the vertices in 
Q are all stored in the vertex normal form, this can be done in constant 
time. 
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Algorithm 6.10. (The reduction algorithm) 

Input: v' G V(T) and Q the output of Algorithm 1 6. 5\ with initial vertex v. 

Output: A tuple (w,j) G V(£?) x T with v' = jw. 

Algorithm: 

1. Let 7o : (v ' = v m , v m -i, . . . , v) be the geodesic from v' to v. Let v.- t be 
the vertex of 7o H closest to t>'. Let r = m — i, this is the distance 
from v' to Q. 

2. If r = 0, we have v' G Then return (v', 1). 

3. If r > 0, we distinguish two cases: 

(a) If Vi is projectively unstable, by a for-loop through the elements 
7 in Stabr(fj), find an element 7 G T such that 71^+1 is a vertex 
of ^. Replace v' by 71;' and apply the algorithm recursively to get 
some pair (w,j) in V(^) x T. Return (w,jj). 

(b) If Vi is projectively stable, run a for-loop through the vertices v 
in Q adjacent to t>j to find the unique v such that either: (i), the 
edge label of the edge from v to V{ is of the form (e, 7) for some 
7 G T with yt(e) = Vi and 70(e) = fj+i, or (ii), the edge label 
from Vi to v is of the form (e, 7) for some 7 G T with o(e) = 
and t(e) = v i+ \. In case (i), replace v' by 7~ 1 f / and apply the 
algorithm recursively to get some pair (w,j) for r y~ l v' . Return 
(w,77 _1 ). In case (ii), replace v' by jv' and apply the algorithm 
recursively to get some pair (w,^) for jv'. Return (w, 77). 

Proposition 6.11. Letv' in T and let Q be the output of Alaorithm \6.5\ under 
the hypothesis of Theorem \ 6.S\ with initial vertex v. Then Algorithm \ 6.1(A 
computes a V -equivalent vertex w of v' and an element 7 G T with 71/ = w. 
It requires 0(n 3 deg(r) 2 ) additions and multiplications in ¥ q where n is the 
distance of v' to Q. 

Proof. In both cases of the algorithm we find an edge label that moves v' 
closer to Q. Since each step of the algorithm decreases the distance d(v',Q), 
the algorithm terminates after at most n steps. From Corollary 19.51 and 
Proposition 14.141 it follows that at step j one multiplies a matrix of height 
(J — 1)| deg(r) with one of height | deg(r). Further one has to compute the 
vertex normal form of a matrix of height at most (j + 1)| deg(r). This takes 

at most (87 + 8j 2 )| 2 deg(r) 2 operations in ¥ q . Summing over j, the asserted 
bound follows. □ 
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Example 6.12. In Figure [2] we give an example of the Algorithm 16.5} where 
q = 5 andr = T(T+ l)(T + 2)(T + 3). We start with as the initial 

vertex v. The adjacent vertices correspond to the matrix , which is 

(\IT 2 ct\lT\ 
q ^ J with a 6 F5. Using 

fl/T 2 0\ 

the algorithm described in Section [7| we compute that I ' ^ j is the only 
projectively unstable vertex and 

'1/T 2 1/T\ fl/T 2 41/TV 



# Homr 











This finishes Step 1 of the algorithm, as depicted in Figure [2J In Step 2 we 
then continue with the eight indicated vertices of level 3. In this case, the 
algorithm terminates after 3 steps. 




projectively stable 
o projectively unstable 



Step 3 

Figure 2: Example: q = 5, r = T(T + 1)(T + 2)(T + 3) 
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Example 6.13. Consider K = F 5 (T) and the two discriminants r\ = (T 2 + 
T + 1) • T ■ (T + 1) • (T + 2) and r 2 = (T 2 + 2) • T ■ (T + 1) • (T + 2). 
Let Ti be the group of units of a maximal order of a quaternion algebra of 
discriminant for i G {1,2}. Then Ti\T has 14 cycles of length 2, while 
r 2 \T has 10 cycles of length 2. Hence these two graphs are not isomorphic. 
This answers a question of Papikian who asked for an example in which the 
lists of degrees of the factors of r and r' are the same but where the graphs 
are non-isomorphic. This is similar to [GN] Rem 2.22] where congruence 
subgroups r (n) and r (n') of GL 2 (A) are considered. 

7 Computing Homr(t>, w) 

Let r, a, A C D = (jf) , t be as at the end of Section HJ recall also that 
ti = 1/T is a uniformizer of K^. 

Lemma 7.1. To compute \fa in = F 9 ((7r)) to n digits of accuracy one 
requires 0(n 3 ) additions and multiplications in W q . 

Proof. Let m = deg(a). It suffices to compute the square root u of the 1-unit 
7i m a to n digits accuracy. This can be done by the Newton iteration in n 
steps starting with the approximation uq — 1. The fc-th approximation is 

Uk = Uk_i — . From the right hand expression one only needs to 

u k — 1 

compute — vr m a; which requires n 2 operations in ¥ q . The k-th digit past 
the decimal point divided by 2 has then to be subtracted from itfe_i. □ 

To state the following result, we define a (logarithmic) height || || on ele- 
ments of A. Its definition will be in terms of our standard A-basis of A, and 
it will depend on this choice. For (Ai, A 2 , A 3 , A 4 ) G A A we define 

st ~\~ %j 

||Ai • 1 + A 2 • % + A 3 • j + A 4 • II := max deg(AA (7.1) 

a ie{l, 2,3,4} 

We also define, as an abbreviation, applied to a matrix or a vector of 
elements in to be the minimum of all the -valuations of all entries. 

Theorem 7.2. Suppose v,v' G V(T) have distance n from Vq = [L(0, 0)]. 

1. There is an algorithm that computes Honir(t>,i/) in time 0(n A ) field 
operations over W q . 

2. All 7 G Hom r (t>,t>') satisfy \\^\\ < n + deg(a)/2. 
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Proof. If v = [L(l,g)] has distance n from Vq, then either 



I = n and deg ; (g) lies in {0, ... , n} or 

n + 1 

I G {— n, —n + 2, — n + 4, . . . , n — 2} and deg ; (g) = — - — , 

see Figure [1] and Remark 13.51 Moreover the path from [L(0, 0)] to [L(l,g)] 
is via £(^p, 0) if I < n and via L(n — deg^g), 0) if I — n. Set n\ := deg ; (g) 
and n 2 := n — ni if / = n and ri2 = n\ — n if Z < n. In Figure [T], the integers 
n\ and n 2 G Z are the coordinates of v from the baseline toward it and along 
the baseline, respectively. Moreover / = n\ + n 2 and g G n^^Ooo = 7r ri2 (9 00 . 
Similarly we define n' x and n' 2 for v' = [L(l',g')} which is also of distance n 
from v = [£(0,0)]. 

/vr' q\ fir 1 ' q'\ 

Let now 7 = I and 7 — I be the matrices in vertex 



normal form representing v and v' respectively. By definition of Homr we 
have 

HomrM') = 7 / GL 2 (0 00 )^7" 1 n T. 

Because u 0O (det(7)) = I, f 00 (det((7 / ) _1 ) = f 00 (det(o")) = for all a G 
GL 2 (C DO ) and f(detT) = {0}, we see that 

HomrM') = n^SGUiO^h- 1 n T, 

where we simply write ^ l ~ 1 )/ 2 for the scalar matrix n^~ 1 '^ 2 ■ 1 2 . By taking 
determinants on both sides and using the fact that n A = ¥ q , we finally 
obtain 

Rom r (v,v') U {0} = ^ l - V)/2 iM 2 {0 00 ) 1 ' 1 n A. (7.2) 

Equation (17.21) can be interpreted in the following way: The intersection 
in the previous line is up to change by conjugation the same as M 2 (Ooo) fl 
n^'- l) / 2 y- 1 A'j. Here M 2 (0 00 ) is the unit ball in M 2 {K 00 ), a i^-vector space 
of dimension 4 and ir^ 1 ~ l ^ 2 y~ 1 A , y is a discrete A-lattice (of rank 4) in this 
vector space. I.e., we need to compute the shortest non-zero vectors of the 
lattice 7P J ~ l ^ 2, y'~ 1 A'j with respect to the norm given by M 2 {0 00 ). If these 
vectors have norm at most one, they form Homr(f,f')- If their norm is 
larger than one, then Homr(f,f') is empty. In particular, the problem can 
in principle be solved by the function field version of the LLL algorithm. 

However, the implemented versions of the LLL algorithm [Hel IPauj need 
an a priori knowledge of the precision by which a has to be computed as an 
element in F 9 ((7r)). This in turn makes it necessary to find a bound on the 
height of the elements in Homr(f,f'), if described as a linear combination 
in terms of our standard A-basis for A. Moreover, |He| IPauj do not give 
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a complexity analysis for their algorithms. To derive these quantities, i.e. 
precision, height and complexity, we proceed as follows. Set 
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Observe that v^gr*-*) > n 2 - ^™ = > -f and that -|/| > -n. 

This implies that v oc (B) > —n. Similarly, using deg(e) < deg(a) and 
computing C~ x explicitly, one finds t>oo(C x ) > — m where we abbreviate 

We now flatten 2 x 2-matrices to column vectors of length 4. Taking the 
explicit form of the A-basis of A from Lemma [4. 171 into account, as well as 
the explicit forms of 7 and 7', the solutions to (17. 2p are the solution of the 
linear system of equations 

CX = Bx, (7.3) 

where A denotes a (column) vector in A 4 and x a (column) vector in O 4 ^. 
The equivalent form A = C~ 1 Bx and the above estimates on the valuations 
of C _1 and B now immediately imply foo(A) > — (n + m). In other words, 
the components of A are polynomials and ||A|| < n + m. This proves (b). 

Next, consider (17.31) in the form B~ l C\ = x. Again by explicit com- 
putation, we have v OQ (B^ 1 ) > — n and v^C) > — max{deg(r), m} =: — d. 
Writing B~ 1 C = Yl'k=-d-^-k'^ k as a power series with 6 M 4 (F ? ) and using 
the bound from (b), equation (I7.3P is equivalent to 



-(n+m) 



X k 7T- k )\ = 



(mod Ot) 



We also expand A = Y^kt^ ^k 7T ^ k as a polynomial in tt^ 1 with X k G ¥ q 4 and 
let Xk and X k be zero outside the range of indices k indicated above. Then 
(17. 3p becomes equivalent to the system of linear equations 

(EST X h - k X k ) =0, h = 0,...,2n + d + m 

in the indeterminates A fc and with coefficients in ¥ q . (Each equation has 4 
linear components.) On the one hand, this shows that we need to compute 
a to accuracy n' = 2n + d + m + 1. On the other hand, we see that using 
Gauss elimination one can solve for the unknowns in 0(n' 2 ) steps where 
each step consists of (An') 2 additions and (An 1 ) 2 multiplications in the field 
W q . Regarding deg(r) as a structural constant and applying Proposition I4.14[ 
the complexity is thus 0(n 4 ). □ 



24 



Remark 7.3. We have chosen vq as a reference vertex in Theorem 17.21 for 
simplicity. Since GL 2 (i^ 00 ) acts transitively on T, one could work with any 
reference vertex. Also, if one chooses vq as the mid point of the geodesic from 
v to v', one sees that the complexity of an algorithm to compute Homr(t>, v') 
is 0{d i ) where d = d(v,v'). Note also that only vertices that are an even 
distance apart can have non-trivial Hom r (t>, v'), because d(v,jv) is even for 
all 7 G T and v G T. 

Remark 7.4. Our implementation of algorithm of Theorem 17.21 uses the 
Gauss algorithm and not LLL. The linear system that needs to be solved 
has An' equations in An + 2 deg(a) variables with n' as in the above proof. 
In practice, deg(a) < deg(r), compare Proposition 14.141 As we shall see in 
Proposition I9.4[ see also Remark l9.6[ we have n < 2 deg(r) — 2 and typically 
< 2 deg(r) — 4. Therefore we have about An' < 22 deg(r) equations in about 
10deg(r) variables. Since the number of vertices of the quotient graph is 
essentially q de ^ r '~ 3 (and q > 3), already deg(r) = 10 is a large value. Over 
finite fields, systems of the size just described can be solved rather rapidly. 

Remark 7.5. To adapt the algorithm of Theorem 17.21 to the generality pro- 
posed in Remark 16.91 at this point requires substantial new code for function 
fields. Using the notation from there, the rings A tend not to be UFD's 
and thus have a more sophisticated arithmetic. Moreover we do not expect 
that one should be able to give explicit simple formulas that describe the 
quaternion algebra D and even less so a maximal order A in it. One could 
work with non- maximal orders A but the quotient graph A*\7~ has typically 
much larger size than A*\T. If one has a reasonably simple description of A 
then an algorithm as in Theorem 17.21 should be doable (certainly in the case 
where S consists of one place only). Also, as far as we are aware of, an LLL 
algorithm in this generality is not implemented. Because of all these still 
open problems, it seems reasonable to present the algorithm here for F 9 [T] 
only. 

8 Presentations of T and the word problem 

From a fundamental domain for the action of T on T together with a side 
pairing one obtains a presentation of T as an abstract group. This has been 
explained in [Sell Chapter 1.4] interpreting T as the amalgam of the stabiliz- 
ers of the vertices of F\T along the stabilizers of the edges connecting them. 
Compare also \P&1\ Thm. 5.7]. 
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Lemma 8.1 ( |Sel[ 1.4.1, Lem. 4]). Let G be a group acting on a connected 
graph X and y a fundamental domain for the action of G on X with an edge 
pairing PE. Then G is generated by 

{g e G e G PE} U {Stah G (v) | v G V(«S)}. 

The relations among the generators of the previous lemma are given by 
\Sel\ § 1-5, Thm. 13] and based on the construction of the fundamental group 
Ti(T,y,S) in |Selt p. 42]. For the group T considered here, all non-terminal 
vertices v of S have stabilizer F* which lies in the center of T. The results 
just quoted therefore considerably simplify and yield: 

Proposition 8.2. Let {y,S, (<7 e )eePE) be a fundamental domain with an edge 
pairing for (T, T) as provided by Algorithm \ 6.5l For each terminal vertex 
v G V(5), let g v be a generator of Stabr(w). Then V is isomorphic to the 
group generated by 

{do} U {g v | v terminal in V(S)} U {g e the edge-label \ e G PE} 

subject to the relations 

g^ 1 = 1, gl +l = g for all terminal v, [g e , go] = 1 for all e G PE . 

In particular go lies in the center of T, as it should. 

Example 8.3. In Example 16.121 the group T is generated by 

{go, Qvii ■ ■ ■ j gv 8 ,gi, • • • > g$} 

with relations 

go = l i9 Q Vi = go, [go,9i] = 1. 

The word problem with respect to this set of generators was already solved 
by the reduction Algorithm 16. 10} compare \Vo\ Remark 4.6]. 

9 Complexity analysis and degree bounds 

In this section we will analyze the complexity of Algorithm 16.51 and obtain 
some bounds on the size of generators of T. We start by bounding the 
diameter of the graph T\T. The idea of using the Ramanujan property 
to obtain complexity bounds was inspired by \KV\ Conj. 6.6]. A standard 
reference is |Luj . 
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Definition 9.1. A A;-regular connected graph Q is called a Ramanujan graph 
if for every eigenvalue A of the adjacency matrix of Q either A = ±k or 
|A| < 2 V / F^I. 

Proposition 9.2 HLu, Prop 7.3.11]). Let Q be a k-regular Ramanujan graph 
on n > 3 verticesu Then 

diam(£) < log^An 2 ). 

Let 



one(i?) : = 



1 if some place in R has degree one, 

q(q — 1) otherwise. 



Lemma 9.3. There is a covering of Q := T\T by a q + 1-regular Ramanujan 
graph Q with 

^ r/ xs 2one(i?) -i-r . 
# V ^ = 7^Vllfe- 1 )- 

Proof. Recall the definitions and formulas for V\ and V q+ i from Theorem 15.61 
If one(i?) = 1, we can choose a degree 1 place p G R. If not we choose 
an arbitray degree 1 prime po- Let r(p ) be the full level po congruence 
subgroup in T. By jTSVl Thm. 1.2] we know that ^ := (m r(p ))\T is 
a Ramanujan graph. Observe that (T fl r(p ))\r = F* 2 if po G R, which 
has cardinality q 2 — 1, and (T fl r(p ))\r = GL 2 (F 9 ) otherwise, which has 
cardinality one(i?)(g 2 — 1). By analyzing the growth of the stabilizers from 
T n r(po) to r, we observe that 

' #V{G) = V 1 + {q + l)V q+1 



one(R) 



peR 

□ 



Proposition 9.4. Suppose V(T\T) > 3. Then 

diam(r\T) < 2deg(r) + 2(21og g (2) + 1 - log q (q - 1)). 



^^The proof in |Luj requires at least one eigenvalue A of the adjacency matrix with 
|A| < 2y/k — 1 and hence n > 3. Also, the assertion is obviously wrong for n = 2 and k 
large. 
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Proof. Let Q = F\T and Q' be the covering from Lemma [9.31 Then 



diam(£) < diam(S') < 21og ? (# V(0')) + log,(4) 

E2 / 2q T-r \ 

< 2io gg (-^ T n^- i ))+ i °g 9 ( 4 ) 

^ peii 

< 41og g (2) + 21og g (-^ I )+21og 9 (n*) 
= 2(21og,(2) + 1 - log g (g - 1)) + 2deg(r). 

□ 

Corollary 9.5. With \\ \\ as in the group T is generated by the set 

{7 e T | || 7 || < deg(a)/2 + 2 deg(r) + 2(2 log, (2) + 1 - log 9 (g - 1))}. 

Proof. By Proposition l8.2[ the group T is generated by the vertex and edge la- 
bels of the quotient graph from Algorithm l6.5l By Proposition l7.2l these labels 
g t have norm \\g t \\ < deg(a)/2 + n, where n is the distance in F\T between 
the initial vertex and the labeled vertex. In particular, n < diam(r\T). □ 

Remark 9.6. If one(i?) = 1, we can obviously subtract 2 + 21og (J (g — 1) 
from the diameter in Proposition 19.41 and subsequently from the bounds in 
Corollary 19.51 In the other case we expect this to be possible as well. This 
should follow by replacing r(p ) by 




(mod po) in GL 2 (F 9 )}. 



Unfortunately we could not find this analog of |LSV} Thm. 1.2] for a congru- 
ence subgroup other than T(p) in the literature although it seems likely to 
hold. 

If this was indeed true, we would obtain the improved bound 

diam(r\T) < 2 deg(r) + 4 log, (2) - 4 log 9 (g - 1). 

For q > 19 it gives diam(r\7~) < 2deg(r) — 4. The nice feature of this 
last bound is that it was assumed in many concrete examples that we have 
computed. 

Proposition 9.7. Algorithm \6. 51 computes the quotient graph T\T in time 

£>((# V(r\r)) 2 diam(r\T) 5 ) l ^ 1 C(g 2dog{r) - 6 ■ deg(r) 5 ) 
in terms of operations over ¥ q . 
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Proof. According to Prop 17.2} comparing two vertices in the algorithm can 
be done in time 0(n 4 ), where n is always less or equal then diam(T\T). 
The list of vertices in each step of the algorithm is always shorter than the 
cardinality of V(r\T), so in each step the number of comparisons is bounded 
by V(r\7~)) 2 . The number of steps is bounded by diam(T\T) and the 
result follows. □ 
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